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A SHORT NOTE ON KRONECKER SQUARE ROOTS 


YORICK HARDY 


Abstract. The results of [1. Ojeda, Amer. Math. Monthly, 122, pp 60-64] 
provides a characterization of Kronecker square roots of matrices in terms of 
the symmetry and rank of the block vec matrix (rearrangement matrix). In 
this short note we reformulate the characterization in terms of rank only by 
considering an alternative to the block vec matrix, provided that the charac¬ 
teristic of the underlying field is not equal to 2. 


Let ® denote the Kronecker product and A the fc-th Kronecker power of a 
matrix A. An m x n matrix A is said to be a fc-th Kronecker root of an x 
matrix M ii M = A. 

Ojeda introduced the notion of the block vec matrix [I] to characterize Kronecker 
square roots and to describe a simple procedure to compute an Kronecker square 
root for real and complex matrices. In particular, the rearrangement matrix in [5] 
plays a central role. Let A be an m x n matrix and i? be an s x t matrix. The 
rearrangement operator Rmxn [2] is defined by 

RmxniA (g) B) = vec(A) vec(B)'^, 

and linear extension. A characterization of Kronecker square roots is given in [T] : 

Theorem 1 ([TJ Corollary 1]). If M is a non-zero m? x matrix and A is an 
m X n matrix then 

(1) M = A® A if and only if Rmxn{AI) = vec(A) vec(A)’^, 

(2) if M = A® A, then Rmxn{AI) is symmetric and has rank one. 

This short note will show that (2) may be reformulated without reference to 
symmetry (Theorem [S]) , provided that the underlying field is not of characteristic 
2. The reformulation of (1) and Theorem [5] follows trivially, as described at the end 
of this note. 

Theorem 2 ([Tl Theorem 2]). If M is an x real or complex matrix such 
that Rmxn{AI) is symmetric and has rank one, then 

(1) there exists an m x n complex matrix A such that M = A® A, 

(2) if A and B are m x n matrices satisfying M = A®A = B®B, then 
B = iA, 

(3) if M is real then there exists an mx n real matrix A sueh that M = A® A 
if and only if tiiRmxniM)) > 0. 

In the proof of part (I) and part (3), Ojeda provides a straightforward procedure 
to determine Kronecker square roots of matrices over the real or complex numbers. 
An analogous method holds using the reformulation of Theorem[T] First, we define 
the “rearrangement” operator R^xm which plays a similar role to Rmxn- 
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Definition. Let Ai, Ak he m x n matrices and j € {1,..., fc}. Define the j-ih 
rearrangement operator 

i?Hxn("4i 0 • • • ® Ak) = vec(Alj) vec(2li 0 • • • 0 Aj_i ® A^+i ig) • • • (g) Ak)'^ 
and linear extension. 

It follows that for fc = 2 we have Rmxn{M) = The 

j-th rearrangement operator is bijective. Consequently M = A if and 

only if 

RlnxniM) = vec{A) vec 

for all j S {1,..., k}. In this case has rank one. 

Remark. The utility of rI^.^„{M) is that it rearranges the entries of M in a 
configuration that is suitable for matrix rank analysis. Other rearrangements or 
unfoldings are equally applicable for this purpose. 

Definition. Let Ai, Ak he m x n matrices. We define R^xn by 

k 

Rl^xniM ®---®Ak)=Y, R^^xu{Al ®---®Ak) 
i=i 

and linear extension. 

Clearly i?^xn i® not injective. 

Lemma 1. Suppose that the characteristic of the underlying field does not divide 
k. The restriction R^^n ■ ^ Rrnxni^) i® bijective, where 

X = I a (g)^ A : where A is an m x n matrix and a a scalar | . 

Proof. Let a (g)^ A G X and (3 (g)^ B G X, then 

if and only if 

kR^L («(2)' a) = kR^^l^ (/3 b) 

and since k ^ Q and R^xn i® bijective we have A = (3^^ B. Thus the 
restricted map R^xn i® injective. The restricted map R^xn i® surjective since 

□ 

Remark. The I^ x 2^ matrix M = (l —1 1 0 0 0 0 O) is not a Kronecker 

cube, but i?fx 2 (-^) i® n rank 1 matrix. Of course, M ^ X in this example. 

With the above definitions and lemma, we are ready to extend Theorem [T] for 
higher order Kronecker roots. 

Theorem 3. Suppose that the characteristic of the underlying field does not divide 
k. If M is a non-zero x matrix and A is an m x n matrix such that M = 
(g)^ T, then R^^^{M) has rank one. 
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Proof. Since for M — A we have that 

RfnxniM) = kvec{A) vec a) , 

and since the characteristic of the underlying field does not divide fc, fc 7 ^ 0 and it 
follows that has rank one. □ 

The connection between Theorem [T] part (2) and Theorem |3] is as follows. 

Theorem 4. Let k = 2 and let M be a non-zero x matrix over a field with 
characteristic not equal to 2. Then Rmxn{M) is symmetric and has rank one if 
and only if has rank one. 

Proof Recall that, for k = 2, R^^^.^{M) = = Rraxn{M). If Rmxn{M) 

is symmetric and has rank one, then 

= 2R™xn(M) 

which has rank one. For the converse, suppose R^^„{M) has rank one. Let 

p 

M = ^Aj ® Bj 
i=i 

be a tensor decomposition of M, where the matrices Aj and Bj are mxn matrices. 
Then 

p 

■Rmxn(^) = [vec(Aj) vec(Bj)^ + vec(Rj) vec(Aj)'^] . 
t=i 

Since R^^„{M) has rank one, we find that the matrices Aj and Bj are all scalar 
multiples of the same non-zero matrix A^,. Thus 

p 

Rmxn{M) = ^ vec(Aj) vec(i?j)^ = cvec(A*) vec(A*)^, 

for some scalar c. Clearly c 7 ^ 0, since if RmxniM) is zero, then M is also zero 
and consequently R^^„{M) is also zero (i.e. does not have rank 1). Consequently 
Rmxn{M) = cvec(A*) vec(A*)^ is symmetric and has rank one. □ 

Theorem [1] part (1) and Theorem [5] parts (1) and (3) and their proofs in [I] may 
be trivially reformulated using Rmxn{M) = when the underlying field 

does not have characteristic equal to 2 . 
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